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Many types of engineering equipment require the flow of 
fluids through screens. Typical applications include filters# 
thickeners, paper*-making machines, and coalescers. Recent studies 
have shown that fine-me^h screens also make an excellent system 
for controlling liquids in spacecraft propellant tanks during 
periods of reduced gravity. Surface tension forces between the 
screen wires and the liquid selectively exclude vapor, thereby 
containing liquid at a desired location. Typical applications 
for these devices are the providing of vapor-free liquid for re- 
starting an engine and for liquid transfer from one storage tank 
to another. To be successful, the devices must be capable of 
being refilled after use. To date, considerable effort has been 
expended in studying basic characteristics such as wicking and 
flow-through pressure drop, but relatively little effort has been 
spent on fundamental understanding of liquid behavior during re- 
filling . 

Several investigators have developed correlations relating 
the velocity V through a screen and the pressure drop AP across 
the screen. Most correlations can be put in a form equivalent to 
the one developed by Armour and Cannon (ref. 1); 



One exception to this statement is the model developed in refer- 
ence 5, for square-mesh screens, which relates AP and V linearly. 
The geometry of the screen weave is represented by: E, the void 

volume ratio; a, the ratio of unit surface area to unit volume; 
b, the screen thickness; D, the pore diameter or particle-reten- 
tion rating; and Q, a tortuosity factor, the ratio of a typical 
fluid path length to the screen thickness. Fluid properties are 
represented by viscosity w and density p. Both a and 0 are 
empirical factors which represent, respectively, viscous and 
inertial contributions to AP . Armour and Cannon's tests gave 
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universal factors for a of 8.61 and 6 of 0.52. More recent tests 
(refs. 2 and 3, for example) indicate that these parameters depend 
upon the screen weave. 

All correlations of the form of equation (1) have been de- 
veloped for uniform, unidirectional’ flow confined by the walls 
of a tube, across a perpendicularly-oriented screen. For space 
applications, it is unli)cely that the refilling liquid will be 
confined; some sort of free-jet or free-flow process is more rea- 
sonable. Tests of the impingement of a circular jet of liquid 
upon various screens under weightless conditions (ref. 4) have 
shown in fact that equation (1) overpx :>dicts the through-flow 
by a considerable margin, a conclusion arrived at by interpreting 
Ap as the average pressure (i.e., the impact force of the jet 
divided by the impact area) on the upstream screen surface. Much 
of the liquid was simply deflected as if by a baffle. 

The analysis presented in this report was undertaken to help 
gain an understanding of the impingement of a liquid jet on a 
screen. Specific objectives were: 

• Develop a mathematical model for predicting the 
through-flow of a circular jet of liquid impinging 
normally or at an angle (45®) onto the surface of a 
fine-mesh screen. 

• Compare and rationalize the model with existing NASA 
data . 

• Define additional experiments to provide an increased 
degree of confidence in the model. 

Computer or wholly nnmdrical models have been excluded, to 
emphasize physical understanding and to enable the model to be 
used economically for design and parametric studies. 




MODELING OF SCREEN THROUGH-FLOW 


Modeling Approach 

Analysis of flow through a screen requires knowledge of the 
driving pressure distribution at all points on the screen. Thus, 
the first modeling task is to derive expressions for the pressures; 
potential flow theory is used for this. Next, a physical model of 
the flow within and through the screen is proposed, using equa- 
tion (1) as a starting point; the formulation of this flow model 
constitutes the main original effort of the analysis. 

Impingement Pressures 

Ideal, irro tat ional , incompressible flow is used to compute 
the pressures exerted on the screen by the impinging circular 
jet, and the screen is idealized as a rigid impermeable plane. 

Most liquids used for space applications have small viscosities, 
so viscous effects are important only in the boundary layer. 
(Boundary layer effects are incorporated in the screen through- 
flow model.) Since the emerging jet has a velocity which is only 
a small fraction of the impingement velocity (ref. 4), its momen- 
tum is negligible in the first approximation, and thus the assump- 
tion of an impermeable screen is reasonable for computing impact 
pressures . 

A previous theory (ref. 6) predicts both pressure and 
tangential velocity distributions for normal impingement of a 
circular jet on an impermeable plane. To obtain predictions that 
are more directly usable, however, an independent numerics] 
computation was made (Appendix A) . These results show that the 
pressure and tangential (radial) velocity distributions can be 
closely approximated by the simple expressions: 
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C = 1 - 0. 221 (r/R . ) 2. 28 
P 3 


( 2 ) 

(3) 
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( 4 ) 


V * 0. 515V . (r/R. ) 
r j ' j 

In developing the through-flow model, it is more convenient to 
approximate these relations by ones that can be analytically 
integrated in certain axpressions. The additional approxima- 
tions are: 


C » 1 - 0.25(r/R.) 
P 1 


V = 0. 5V . (r/R. ) 
r 3 3 


(5) 

(6) 


There is only a small difference between equations (3) and (5), 

2 2 

or between (4) and (6). Furthermore, the impact force, TTpv^R^, 

2 

and the relation between C and V., namely C » 1 - (V /V.) , are 

P 3 P r 3 

given exactly by both approximations. Incidentally, positive 

AP implies that liquid will flow through the screen. Thus, 

equation (5) indirectly verifies the conclusion of reference 4 

that the emerging jet has a diameter about twice that of the 

impinging jet, since C >0 for r <• 2R. . 

P 3 

Previously developed models for oblique impingement apply 
only to two-dimensional or rectangular jets (reJs. 6 and 7). 
Numerical solutions, even for one angle of obliquity, are im- 
practical. Consequently, an approximate procedure described 
in Appendix A is used to predict the pressure distribution. To 
get integrable expressions for C^, analytical expressions are 


required; these approximations are: 
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1 - 0.25(r/R^) 


-90“ < 0 < 90“ 


(7) 


C = U - (r/R.)^l cos^e + n - 0.25(r/RJ^] sin^9 
P 3 3 

90 * <0 <270 


= 0.5V^ ^ ^ 2 ^3 


(8) 

(9) 
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The coordinate system is shown in the sketch. Equation (7) 
neglects a region of small pressures shown schematically by 
the shaded area (see also discussion in Appendix A and Figure 
A-3) . Also, for the term multiplied by cos^6 in equation (8 ) , 


r must be limited to values less than . The total force, 
2 2 

TTpv R cos6, is overpredicted by equations (7) and (8), but 
by less than fifteen percent. 



Predictions of Unmodified Model 


The applicability of the duct-flow relation between pres- 
sure and velocity, equation (1), can be determined by substituting 
equation (5), for normal impingement, or equations (7) and (8), 
for 45® oblique impingement into equation (1) . The solution for 
V is 



2 2 2 

where A * a(a bQ/e ) and B = B(bQ/DE ). Equation ( 0) gives the 
velocity of the emergent jet as a function of r since Cp depends 


5 



on r. To compare this equation with test results# it is neces** 
sary to make the assumption that the average velocity over the 
area of the emergent jet corresponds to the measured velocity 
with which the jet leaves the screen as a growing geyser (ref. 4). 

The average velocity for normal impingement is derived by 
integration over the range of r for »hich AP ^ 0# namely 
0 < r < 2R.. The result is 



( 11 ) 


This expression overpredicts the test data of ref. 4 by about 

loo percent# on average. Reference 4 proposed an empirical 

correlation for V which essentially amounts to the use of an 

a V 

average C of 1/8, rather than the 1/4 resulting from equation 
P 

(5); that is# the effective pressure is supposed to be one-half 
the actual average pressure. Predictions f^*om this correlation 
are within ±15 percent of the test data (ref. 4). There is no 
apparent underlying reason that can explain a reduction of AP 
by a factor of 2; consequently, the correlation may not be 
reliable for other test data. In any case it does not aid in 
understanding screen-flow mechanics and so cannot be extended 
to other kinds of potential applications involving unconfined 
flows unless first verified by experimental data. 

The duct-flow model also overestimates test data for 45® 
impingement by about 100 percent. 


Impinging Jet Screen-Flow Model 

Since the duct-flow screen model does not give reliable 
predictions of the through-flow for jet impingement, a new or 
modified model is required. It seems likely that liquid flowing 
along the plane of the screen or tangentially within the screen 
is the reason for the failure of the duct-flow model. The duct- 
flow model assumes a more or less straight-through path for the 
liquid, but liquid will be deflected into the plane of the screen 
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for all unconfined flows. There is an additional pressure drop 
caused by the increase in the flow^path length resulting from the 
flow deflection which is not included in the duct-flow model. 

This additional AP can be incorporated into a screen-flow model 
in several ways. The original basis of equation (1)/ namely that 
the pressure drop is caused by a combination of viscous flow around 
the screen wires and inertial losses through parallel intercon- 
necting channels between the wires, could be altered by assuming 
the screen to be a porous medium. Or the basic assumptions of 
equation (1) could be kept, but the trajectory of the flow altered 
to account for the in-plane component. 

For closely woven screens such as the twilled Dutch weave 
shown in the sketch, it may be plausible to idealize a screen 
as a sort of porous medium. It is further plausible that the 



internal resistance for tangential flow in the plane of the screen 
is less than the resistance for flow across the thickness; i^ 
fact, according to the theory developed in Appendix C, this is 
required in order to make the predicted through-flows agree with 
the test data of reference 4. Simple calculations (Appendix C) 
show that the in-plane resistance must be less than 1/1000 of 
the through-flow resistance to obtain correlation between test 
and theory. A flow resistance ratio of 1/1000 corresponds, 
roughly, to average pore diametern in the plane of the screen 
that are at least 30 times larger (for inertial losses) or 1000 
times larger (for viscous lO'?/ .") tV. n the pores in the through- 
flow direction. These values seem too large to be reasonable. 
Thus, the porous medium model is probably not physically realis- 
tic . 
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Models which retain the basic assumptions of equation (1) 
can ;*ccount for the non-direct flow path by increasing the tor- 
tuosity factor. Q. as shown conceptually by the sketches. The 
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modified tortuosity factor can be estimated by assuming that the 
flow paths are inclined at an angle which is determined by both 
the through-flow velocity V and the tangential velocity The 

2 2 1/2 

flow path length is then H * bQ{V + V^) /V. Note that the 

length varies over the impingement since both V and are func- 
tions of r. For this model, equation (1) is modified to account 
for jet impingement dynamics by replacing the flow length bQ by 
2 2 1/2 

i, and V by (V + V^) , Thus, the expression 



( 12 ) 


is the pressure-velocity relation for this model. Although this 
model predicts a reduction in V for a given AP when compared to 
equation (1), it has several shortcomings. First, the flow paths 
cannot be inclined at all angles throughout the range 0® to 
90®. Equation (12) requires this, however, since: = 0 and 

V > 0 near the stagnation point of the jet, which gives an in- 
clination of 0®; and V = 0 and V >0 near the emerging jet 

s 

periphery, which gives an inclination of 90®. Second, for large 
inclination angles, the viscous shear forces at the interface 
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between the screen and the impinging jet have significant com- 
ponents in the flow-path direction; this, as well asAP, helps 
drive the flow. Equation (12) neglects the shear force. It 
is evident that by doing so an unrealistically large AP is 
required near the periphery where the through-flow is small 

(V * 0) , since V is large there. The contribution to the total 
s 

through-flow of the peripheral parts of the jet is small, however, 

and equation (12) may model the important central parts of the 

through-flow fairly accurately. The overall effect of in 

this model is unrealistic, nonetheless, because of the unbounded 

behavior of the terms V^/V and V'/V as V 0. 

s a 


To correct these deficiencies, a further modification is 
proposed. As shown in the s)cetch (for a twilled Dutch weave) 
the shute wires which extend into the boundary layer intercept 
a part of the tangential velocity. The intercepted tangential 
velocity is where 4) is a "blockage" coeffic*ent which may 
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depend on both the screen and the jet geometries. As bclore the 

tangential velocity 4>V , combined with V, effectively selects a 

s 
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flow path inclinad at tha angle arctan (V/(|)Vg). The coefficient 
0 accounts for a number of factors: not all flow path angles 

are possible; the inclination of the path changes as the fluid 
passes through the screen, eventually being normal to the rear 
surface; and the screen weave is not axisymmetric about the 
stagnation point of the jet. The model so far proposed is simi- 
lar to the previous model, except for the factor <}>, it is 

further postulated that the tangential velocity enters, exits, 
re-enters, re-exits, etc., the upper layer of wires as the flow 
is dragged along. Thus, the velocity passing through the 3 creen 

2 221/2 

is only the through-flow V, and not (V + 0 V )' . In this 

s 

model the tangential velocity helps select a flow tube inclina- 
tion but does not otherwise affect the pressure drop. As in the 

2 221/2 

previous model, the flow tube length is it = bO(V + (j) V^) /V. 

The relation between pressure drop and through-flow velocity for 
this model is 



or in terms of the constants A and B 


AP = Ay\/v + (0V ) + Bpv\/v + ($v ) 


(14) 


This model, like the previous one, has the disadvantage of r.ot 
model; ’vj large flow angles very well. In fact, it will be seen 
that equation (13) or (14) only makes mathematical sense over 
a surface area for which tlie r-coordinate is less than some R 

max 

2 2 2 

such that V > (AP/Au) - V The reason for this is acain 
— s 

neglecting the viscous shear. It ought to be recognized, 
nonetheless, that this model is valid over a large part of the 
impingement area. As V ^ 0 near the periphery# the effect of 

3 2 

V enters as <J)V , and not as the unboundea terms V /v and V /V. 
s s s s 



Furthermore, for reasonable values of (less than one), numer- 
ical results to be cited later show that the through-flow in the 
region r > R ought to be negligible because dP ^ 0. 

The boundary layer velocity and thickness are estimated 
from previous results for laminar unidirectional flow along a flat 
plate (ref. 9). All the tests of reference 4 had laminar 
boundary layers. Neglecting the radially diverging character of 
the boundary layer is an approximation, but the correct physical 
parameters are introduced and the empirical parameter <p can also 
account for the difference between radial and unidirectional 
flow. The boundary layer thickness is 5 = 5 /pr/pV^ which reduces 

to 6 = 6.97 /pR . /pv when equation (6) is used for V . Since 
13 r 

the screen penetrates into the boundary layer by about one-half 

a wire diameter (b/6 for a twilled Dutch weave) the effective 

slip velocity V can be shown to be 
s 



(15) 


This expression has been derived by the linearized approximation 
b 

V = — (— T ) for 2 = 0, where is the laminar boundary layer 

s 6 dz BL 

velocity profile (ref. 9). The numerical constant 0.017 is 
retained primarily to get the correct order of magnitude for 4» . 
(Only the quantity ({)V^ enters the modeling expressions so the 
pure number 0.037 could in principle be absorbed into <j> . ) Equa- 
tion (15) has been derived for a twilled Dutch weave, but is 
valid within a factor of two for other weaves; this extra factor 
can be absorbed into 4> . 

Equations (14) and (15), together with the imp ,gement pret- 
sure expressions, form the proposed mathematical model relating 
the emergent velocity to the impinging jet characteristics. 


Solution for Emergent Velocity 

Equation (14) must be averaged over the emerging jet area to 
predict the average through-flow velocity and the flow rate 

through the screen. The integration needed to obtain this aver- 
age cannot be done exactly by analytical techniques. Since numer- 
ical integration first requires a value of (p to be assigned, it 
is inconvenient to use for verifying the model with experimental 
data. Instead, an approximate analytical solution is developed 
in Appendix B. For normal impingement, the integrated result is 
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For oblique impingement at 45®, the average through-flow velocity 
is 
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With respect to exact integration, equation (17) is more approxi- 
mate than equation (16). 
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for normal impinge- 
> 0 when AP > 0) , is 


A similar but : e complicated expression can be derived for 
oblique impingement. 

The degree to which the mathematical model of the through- 
flow represents reality can be inferred from comparison with test 
data. A "good" comparison requires that (1) a single value of <p 
correlates all test data for a given screen weave and impingement 
angle, regardless of fluid properties or jet velocities and diam- 
eters, and (2) R not be too much smaller than 2R . . If these 

max 3 

requirements are met, then it can be reasonably concluded that 
the proposed model does incorporate all the important physics of 
the flow, is valid, and can be used reliabily to predict through- 
flows . 


RESULTS AND DISCUSSION 

A substantial amount of test data is available from experi*> 
ments conducted at NASA-Lewis Research Center. Results for normal 
Impingement have been published (ref. 4) , but oblique jet results are 
available so far only as high-speed 16 mm motion pictures. The tests 
were conducted in a weightless environment, using a drop tower. Very 
little liquid penetrated the liquid at low impingement velocities or 
for very tightly woven screens, and any liquid that did tended either 
to puddle or to form a small non-growing geyser on the back side. 

This type of behavior was correlated with a Weber number criterion, 
which can be used to predict the minimum velocity needed for sub- 
stantial through-flow. At higher velocities, some of the jet flow 
penetrated the screen and formed a continually growing geyser. The 
steady state through-flow measured in these tests constitutes the 
data used to verify the proposed screen-flow model. 

Verification of Model 

Twilled Dutch Screens. NASA test results for normal and oblique 
impingement are shown in Table I. Two other screens (325 x 2300 and 
200 X 1400) were also used, but through-flow was not observed with 
these screens for the tested conditions. Geometrical constants 
for the screens and physical properties of the fluids are given in 
Tables II and III. All impinging jets had diameters of 0.625 cm. 

For a given screen and set of test conditions, equation (16) 
for normal or equation (17) for oblique impingement was used to pre- 
dict the average through-flow velocity as a function of the block- 
age coefficient Typical results for normal impingement are 

shown in Figure 1. As (p increases, the through-flow decreases, 
which is the desired behavior. The range of (J>-values needed to 
correlate the test data is fairly small, regardless of liquid pro- 
perties or jet velocities, thus tending to validate the assumptions 
of the model. (The only exception is the TCTFE test using V^ « 125 
cm/sec; a value of 0 = 1.1 is required to predict the through-flow 
velocity of this test. The correlation proposed in reference 4 
also overpredicts this test point, by about 100 percent. It is 
believed that the data for this test is in error.) There is, at 
least for this screen, a tendency for the best correlation of (P 

1 4 


TABLE I. VALUES OF BLOCKAGE CONSTANT, (p , FOR TWILLED DUTCH SCREENS 
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EMERGENT JET AVERAGE VELOCITY, V, cm/sec 



BLOCKAGE COEFFICIENT, 0 

Figure 1. Emergent Jet Average Velocity as Function of Blockage 
Coefficient. 200 x 600 Screen and Normal Impingement 



considering only the TCTFE tests to be slightly greater than the 
best correlation if only the ethanol tests are considered. One of 
the premises of the model is that should not depend on fluid pro- 
perties* but only on screen and jet geometrical factors. The ob- 
served behavior of is not large enough to invalidate this premise 

Typical results for oblique impingement are shown in Figure 
2. There is somewhat more scatter in the value of (p needed to 
predict each test than there is for normal impingement. At least 
two causes contribute to this. First# the math model itself is 
less exact because of the approximate pressure distribution used 
as its basis and the number of mathematical approximations needed 
to derive equation (17). Second# twilled Dutch screens do not 
have the same geometrical properties in all lateral directions# 
but the tests were not all conducted with the same relative orienta 
tion of the screen weave and inclined-jet axis. In addition to the 
scatter, the TCTFE data seem to be correlated by larger values 
of $ than ethanol, similar to the behavior shown in Figure 1. The 
model fails to predict the test results for ethanol when » 281 
cro/sec. In fact even 0=0, which corresponds to a duct-flow 
model and a, jet-impingement pressure distribution, underpredicts 
the test result slightly. Consequently, this test is neglected in 
determining the best correlating value of 0. An analogous situa- 
tion (i.e.# the best value of 0 being much smaller for one test 
of a group) occurred for the 200 x 600 screens with ethanol and 
* 281 cm/sec. Both of these tests used jet velocities which 
resulted in little through-flow. 

The overall best correlations for 0 are shown in Table 1. 
Predicted values of using these values of 0 are also given, 

except for the two cases mentioned previously for which the 
model seems to fail. Several conclusions are evident from the 
tabulated results. Blockage constants for a given screen are 
always larger for normal Impingement than for oblique impingement. 
Further# the correlation between the measured and predicted 
emergent-jet velocity is better for normal impingement than for 
oblique impingement; this is the result of the smaller variation 
in the best 0 from one test to another for normal impingement . 


EMERGENT JET AVERAGE VELOCITY, V, cm/sec 
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VELOCITY 

o ACTUAL DATA 


0.2 0.3 

BLOCKAGE COEFFICIENT, 0 

Figure 2. Emergent Jet Average Velocity as Function of Blockage 
Coefficient. 165 x 800 Screen, 45* Impingement 




The largest value of the parameter which deternlnes 

R , is 0.52 » a value that occurs for the 80 x 700 screen and 
uiax 

the smallest ethanol jet velocity. From equation (16), R » 

max 

1.25 R. for <t>iT ■ 0.52. For the bulk of the test data, is 

j 3 3 

smaller than 0.52 (a value of 0.35 is more appropriate), so the 

model is mathematically valid over most of the impingement area 

for most of the tests. 

The overall quality of the correlation of the model for 
normal impingement is shown in Figures 3 and 4, for the 200 x 600 
and 80 X 700 screens respectively. Comparisons of predicted and 
measured velocities are about as good as that shown in reference 
1 as the basis for the duct-flow model. 

Plain Square Screens . Results from two normal-impingement 
tests using a 400 x 400 plain square weave screen were also 
compared to the model. The predicted values for 0 s o are close 
to the test results, as might be expected because of the open 
weave of these screens. Both tests employed ethanol. For = 
214 cm/sec, the predicted value from equation (16) with (}> « 0 is 
28.4 as compared to the test value of 26.5; and for = 246 
cm/sec, the prediction is 36.3 as compared to the test value of 
30.3. For such relatively open screens, the concept of a 
blockage coefficient has little meaning, and in fact numerical 
results show that the predicted through-flow velocity decreases 
by about only ten percent over a range of 4> from zero to one. 

Note that in making the predictions the general values of 
a « 8.61 and 6 = 0.52 (ref. 1) were used since specific data 
for a 400 x iOO screen is not available. Considering the spread 
of the data r-ihown in reference 1 about this recommended general 
correlation, the agreement between test and model predictions for 
a 0 is very lood. 

Percentage of Liquid Passing Through the Screen 

The total flow rate through the screen can be calculated by 
multiplying the average through-flow velocity by the cross- 


2 0 



IMPINGEMENT VELOCITY, cm/sec 

Figure 3, Comparison of Data and Model, Using Best Correlation for Blockage 

Coefficient, 200 x 600 Screen, Normal Impingement 



ETHANOL DATA 



IMPINGEMENT VELOCITY, cm/sec 

Figure 4. Comparison of Data and Model, Using Best Correlation for Blockage 

Coefficient, 80 x 700 Screen, Normal Impingement 



section area of the emergent jet. For normal impingement# the 
emergent area is four times the impingement jet area, or 4 ttr 2. 
For impingement at 45 degrees# the area can be calculated from 
the equation of the emergent jet periphery, r*R^ (2 + cos6) , as 

2 

explained in Appendix A. This area is 4.5 iTR^ . 

Figure 5, a nondimensional plot of equation (16), can be 

used to obtain rapid estimates of the average through-flow 

velocity. Thi&< plot also shows that the inertial pressure drop 

2 

(e.g.# the term proportional to pV in equation {!)) affects the 
through'flow only for high impingement velocities < 3). For 

larger values of (smaller V^) the through-flow is a function 
solely of viscous forces. (The parameter tt^ represents the iner- 
tial characteristics, B# of the screen.) 


NON-DIMENSIONAL EMERGENT VELOCITY, VIT./Vj 



PARAMETER 

Figure 5. Plot of Equation ( 16 ) 
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CONCLUSIOHS AND RECOMMENDATIONS 


Physical and mathematical models have been proposed for an 
unconfined liquid flow Impinging upon a fine-mesh screen. The 
physical basis of the model is the relation between through-flow 
velocity and pressure drop previously developed for confined duct 
flow across a screen. This relation has been modified to account 
for the increase in flow-path length through the screen which 
results from the deflection of an unconfined flow by the screen. 

A detailed mathematical development was presented for a circular 
jet impinging either normally or obliquely (at 45 degrees) on a 
screen. An important result is that a new parameter, the block- 
age coefficient, is required to model the impingement flow 
accurately. Empirical values of the blockage coefficient were 
determined by comparison with test data from previous studies. 

The following conclusions were drawn: 

1. The model verifies the previously observed experimental 
result that the diameter of the emerging jet is about 
twice as large as the impinging jet. 

2. Blockage coefficients depend only on screen geometry 
and jet impingement angle. Coefficients for normal 
impingement are somewhat larger than for oblique 
impingement . 

3. The proposed model correlates normal impingement results 
more accurately than oblique impingement results. This 
is probably not a basic shortcoming of the model but a 
result of approximations made in the mathematical devel- 
opment and scatter in the oblique-impingement data. 

5. Predicted through-flow velocities for normal impingement, 
using the recommended correlations for the blockage 
coefficient, generally agree well (within ±15 percent) 
with the test data. Predictions for oblique impingement 
are less accurate. 

Additional experiments are recommended to provide an increased 
degree of confidence in the model. 

1. Tests similar to previous NASA studies should be con- 
ducted with a wider range of liquids to investigate 

2 5 


possible dependence of the blockage coefficient on 
liquid properties. 

2. Further oblique impingement tests should be conducted, 
taking care to preserve the same relative orientation 
of jet and screen-weave geometry from test to test. 

Other tests should be conducted with the screen weave 
rotated 90 degrees. 

It would also be desirable to run experiments that directly 
verify the basic premise of the model, namely that the increase 
of the fluid path length is the cause of the failure of the duct- 
flow model. Perhaps such tests could be conducted by placing a 
screen in a duct at an angle to the upstream axis. The duct 
walls should be removed downstream of the screen, to eliminate 
wall effects. Measured pressure drop across the screen for a 
given duct velocity could be compared to tests when the screen 
is inserted perpendicularly to the axis. Any increase in pres- 
sure drop can then be attributed to the increased flow path 
length of the screen inclined to the upstream flow. 
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APPENDIX A. IMPINGEMENT PRESSURES 


Normal Impingement 

The only physical parameters needed to describe the ideal 
flow of a circular jet impinging normally against an impermeable 


surface are: fluid density p, jet velocity far upstream , and 


jet radius far upstream R^ . Nondimensionalizing the ideal flow 


equations by these quantities allows all cases to be given by 


one solution. Velocities are nondimensionalized by V ^ , distance 


by R j , and pressures by 1/2 • In terms of a nondimensional 


stream function the nondimensional velocities are 


Lik 

r3 z 


9 ijj 


r9 r 


(A-1) 


Conversion of mass is automatically satisfied by ip , but the 
condition of irrotationality requires 


_9J) 

.2 r9 r 
d r 



(A-2) 


The boundary conditions that solutions of equation (A-2) must 
be subjected to are 


JJi 
r9 r 

= -1 

r < 1 , z -> + ® 

9ip 

r9 r 

= 0 

o 

A| 

o 

II 

H 

9 ^ 
r9 z 

» 0 

K 

II 

O 

N 

|v 

o 

i> = 

0 

on jet and z = 

i|j = 

-1/2 

on free boundary 


(A-3) 


(A-4) 


(A-5) 


(A-6) 


Equations (A-6) are a consequence of: (1} the total flow enclosed 

2 

by the free boundary is irR^V^; and (2) the flow, in terms of the 
stream function, is -2ir times the difference between the values 
of ip at the rigid boundaries (or jet centerline) and the free 
boundary (ref. 10). The negative sign in equations (A-3) and 
(A-6) implies that the jet approaches with a velocity directed 
downward . 

In setting up the numerical solution, a first guess of the 
free boundary shape was made and a solution for ij/ computed. The 
equations were written in finite difference form, using a non*> 
uniform grid of 195 points concentrated near the free boundary to 
improve resolution. A numerical iteration routine was used which 
does not require matrix inversion, so a desk-top programmable 
calculator was sufficient. The accuracy of the first guess in- 
corporated the fact that the height of the deflected jet along 

2 

the wall approaches R^/2r as r ->• ». Several iterations were 
needed to obtain ±5 percent accuracy for the free-boundary ve- 
locity V.; the corresponding accuracy of the surface pressure 

^ 2 
integral differed by less than 3.5 percent from irpR^V^ . 

Computed results for the tangential velocity V along the 

^ 2 

surface and surface pressure coefficient = (p - p^)/(l/2 pV^) 

are shown in Figure A.l. Also shown is the empirical correlation 

2 28 

C = 1 - 0.221(r/R.) ’ . This correlation fits the computed 

curve very well for r < 1.8 R^ and predicts the impact force 

2 2 , . 

TTpR.V. exactly; for r > 1.8R. the fit is not good, but C is so 
j 3 3 P 

small here that the error caused in through-flow calculations is 

negligible . 


Oblique Impingement at 45® 

Numerical solutions for the oblique impingement of a circular 
jet are beyond the scope of the work, and exact solutions do not 
exist. Results for a plane jet impacting at 45® (ref. 7) are 
shown in Figure A-2 (the curve labeled 2 - D , y “ 45®) along with 
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the axlsymmetrlc case for normal impingement and for a plane jet 
(ref, 8) impinging normally (the curve labeled 2-D, y “ 90®). 

For normal impingement, the axisymmetric pressure distribution 
is "fuller" than for the plane jet; and for plane oblique jets, 

# 

the stagnation point moves rearward by about one jet radius. 

Using these results for guidance, the pressure distribution of 
a circular jet impinging at 45® is hypothesized to be as shown in 
Figure A-3. The curves labeled C„, , and are the pres- 

sure coefficients in the coordinate directions shown in the sketch 
below. The distribution was obtained from the inclined 



plane jet result shown in Figure A-2, by making it fuller by about 
the same amount as the difference between the circular and plane 
jets for normal impingement. The Cp^ distribution transverse to 
the jet axis is identical to the pressure for normal impinge- 
ment of a circular jet. At intermediate locations, the pressure 
is assumed to vary smoothly between Cp^ and Cp^f between Cp^ aJ^d * 
Cp 2 ? the proposed distribution is 
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g at 45 




Cp “ Cpj^ cos^e + Cp^ sin^e 


-90* < 0 < 90' 


(A-7) 


Cp 2 cos^e + Cp 2 sin^6 


90® < e < 2/0' 


(A-8) 


Equations (A-7), with the C 's taken from Figure A-3, predict 

P 

2 

an impact force of 0,69 wov^ , which is very close to the 

2 2 

correct impact force 0.707 Trpv. R ^ . 

The region AP > 0 extends to r 3R^in front of the stagna- 
tion point, to r ^ R^ behind the stagnation point, and to 
r 2R^ on either side. A reasonable assumption for the emergent 
jet area is, therefore, 

cos 6 ^ (A-9) 


r < 2R 


. + I cos e ^ 


In the through-flow model, the difference between and 

Cp 2 is neglected, and equation (A-7) is replaced by 

Cp ' Cp2 -90* 101 90* (A-10) 

This is done primarily because an expression of the form 
2 

1 - E (r/R^) where E is a constant, cannot be made to fit the 
Cpj^ curve with any degree of accuracy. However, a quadratic 
form is needed to perform the analytical integrations described 
in Appendix B. The neglected pressures, in any case, are not 
large. The approximate curve-fits for Cp^ and Cp^ discussed 
in the text overpredict the total force by fifteen percent, as 
mentioned earlier. This error incurred by neglecting Cp^^ is 
masked by this approximation. 
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APPENDIX B. DERIVATION OP THROUGH-PLOW VELOCITY EXPRESSIONS 

Normal Impingement 


For normal impingement, the model developed in the section 
titled Modeling of Screen Plow resulted in an expression relating 
the through-flow velocity V to the pressure AP : 

AM \/v^ + (4>V ) BpV Jv^ + (4)V ) ^ '= AP = ^ pc V^ (B-1) 

y s V s 2 p 3 

with 


C =1-0.25 (r/R. ) ^ 
P 3 


(B-2) 


V = 0 .037V . (r/R . ) Jpb^V./pR. (B-3) 

s 3 3 3 3 

To determine the average through-flow velocity and the total flow 
rate, equation (B-1) must be solved for V and the solution inte- 
grated over the through-flow area. Unfortunately, equation (B-1) 
cannot be solved exactly, except by numerical methods which requir 
a defined value for (p . Therefore, the method of successive approx 
imations, V = + AV + . . . , is used, where the first approxi- 

mation V^ is the solution when the inertial term (the smaller of 
the two for fine-mesh screens) is neglected. The result for 
is 


'^l " y (PCpV ./MA) - (4>Vg)^ 


1 


(B-4) 


Equation (B-4) makes sense only for r < R such that C V./AM 

max p j 

2(pV . From equation (B-2) and (B-3), the limiting r is 
s 


r<R =2|\/l + (4>^o ) *■ - 07T 
- max * “ ^ 
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(B-M 


The linear correction, Av, to the first approximation is 
derived by substituting V ■ + AV into equation (B~l) and 

neglecting higher powers of Av. 



(Note that the denominator of the fraction is always greater than 

zero for r < R ) . 

— max 

By comparing V “ + AV to the exact solution for 4> » 0, it 

has been found that the accuracy of the approximate solution is 
±5 percent. Therefore, additional corrections are neglected, 
and the solution of equation (B-1) is assumed to be 



The average velocity is 


V 

av 



2TTrV dr 


(B-8) 


^ 6 


where the flow-through area it( 2R^) is used to compute Ex- 

tending the upper limit of integration to 2 R^ gains nothing since 
for r > R the model breaks down, and it has been argued pre- 

iiiQ X 

viously that only a negligible amount of flow passes through the 
area r > R 

max 

After inserting equation (B-2) for C and (B-3) for V , the 

P s 

first term on the right of equation (B-7) (i.e., V^) can be 

integrated : 



^ <l + 2(4.lij)^+4{(j)7T^)^(l + (4>1T^)^]£n 


4>tt 




1 + (4*1^3) 


(B-9) 


The second term (i.e., AV) cannot be integrated until a further 
approximation is made. Writing this term in the form 



where x = r/R., it can be compared with the simpler expression 



The numerators of these two expressions are equal. Comparisons 
of the denominators are shown in Figure B-1 for various values 
of ({'113 = Y* The simpler expression is an excellent approxima- 

tion except near = R (where the curves break upward Kbarpl'O . 

max 

It has the additional virture of being integrable. For these 
reasons, equation (B- 6 ) is replaced by the simpler 1 ntegrable 
expression. The integral can be evaluated by algebraic division 
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[l-(X)j] dO 3niVA 
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it 



1 


R /R. 
max 3 


(2ftr^) 


•rr^x +4 it2{1 - it2)x + 


ISiT^x 


4(1 ^-TT^) - IT^X 


dx 


(B-iO) 


or 


/Av\ 11 9 r 9 1 

\v“) = - sir 2(<j>TT3 )^ +ir 2 

' 3 / av If 

r , r , [ti + 2 <flt 3 ('/l7((!,ir )2-,j,,,J 

+ 4(|)TT2 |2 (4)TT3) ^ + itjJ y 1 + ((|)TT3) ^ - 27 t 2 fcn ^ = 


The total average through-flow velocity is, therefore. 


V 


H*1I +2(4iit,)^ 
V. 4 tt \ 2 2 ^3 

3 1 


2 + 2((|)TT^) + tt ^ 


+ 4((j)1T^) 


(B-11) 


"1 <t>Tf, 


1 + 


- 2(j5TT, 


2((j)7T^) tlT^ 


1 + (4>''T3)^+ 


7- 

TT^ + 24)77^ ((/TT {(JtTT^)^ - 4>TT3) 


1 + (<}>tT3) 




(B-12) 


Oblique Impingement at 45“ 

The same general procedure is used for the oblique impinge* 
ment, although the integrations cannot be perfomed so neatly 
because of the angular dependence of C^. 

Over the area -90“ ^ ^ s. r '> 0 , the e ur iiace pressure 

for oblique impingement is assumeu to be the sam > as for normal 
impingement in the approximation used here. Thus, the contr— bu 
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tion to V of the flow in this region is merely one-half of 
av 

equation (B-12) . 

Over the area 90® 5. ® 5. 270® (the area behind the stagnation 
point) the pressure coefficient Cp 2 is defined only for r £ R j , 
and Cp 2 defined only for r £ 2Rj (1 + (1/2 cos 0)]; see equation 
(A-10) . Because of this, the various integrals must be handled 
separately. For example. 


3TT/2 R 

/ / 

tr/2 0 


max 


V^r dr d0 


37T/2 R. 

- / / 

TT/2 0 ^ 


r dr d8 


3TT/2 2R^ [1 + (1/2 cos 0) 1 

+ / / V.r dr d0 

7T/2 Rj 


(3-13) 


In the first integral Cp = Cp^ cor *^-0 + Cp^ sin^6 = 1 - x^(l - 3/4 sin^6) 
2 > 

“ 1 - Cx , where ^ * 1 - 3/4 sin“0. This is of the same form as 
the expression used for for notraal impingement. Thus, after 
one integration the first ntegral is 




tt/2 




(1-0^- 


+ 2 I 


+ 


_1 

2C 




Jin 


(<|)Tf^)^ 


_((*)t^) -2^(C~1) -2C ya-^)“- (#3)^] 


n lci0 


(B-14) 


^ — 3/2 

In this expression, ff = (1+1/2 cos 0) ' ] ; the angular 

dependence is a result t the tangential velocity relation V = 
0.037 (r/R^)[(pb Vj)/WKj] ' (1 + 1/2 cos 6) * , which in turn 

result s from the assumption that the periphery of the emergent jet 
is given by r ^ 2Rj(l + 1/2 cos 0) as discussed in Appendix A. 
Since (1 - = ( 3/4 ) sin^O , ti.a integral must be evaluated 

numerically in two parts. The first part covers the range 
(1 - C) >> (say, sin^6 > 0.1); the second covers the range 


4 0 


sin^6 ^ 0.1, for 
first range, the 

'v 2 

powers of ( 4 *^ 3 ) 


which (1 - C) ^ “ 
term [ (1 - 3) ^ - 


to give 


(1 - K) 


1 


2 - 

4> ir^ is close to zero. In the 
(j,2^2]l/2 jjg expanded in 


4>^3 


kth- 

8 1 - C 


Powers of greater than the second are discarded. With t>^ese 

2 

assumptions, expression (B-14) for sin 6 > 0.1 reduces to 



(B-15) 


This expression can be evaluated numerically since <p enters only 
as a multiplicative constant. 

For sin^S j< 0.1, the remaining part of the first integral is 
in ^/o^ 1 

||-+ j (4>7F^)^U - cos 9)^ (l-|-sin^6) ^ | d6 (B-16) 



For convenience the limits of integration in both expressions 
(B-15) and (B-16) have been taken in the quadrant 0 to ir/2 , and 
the term 1 + cos 0/2 has been replaced by 1 - cos 0/2 to account 
for the fact that cos 6 is always negative in the actual range, 
it/2 to 31 T/ 2 . 

The average through-flow velocity contribution of the first 
integral of equation (B-13) is obtained by dividing the integral 
by 7 Tr 2/2. Thus, 



37T/2 R. 

ff 


7T/2 0 


V. 

V.r dr d0 = 

1 


1 2 
i - 0.123{(j5ir2) 


4 1 


(B-17) 


The second integral of equation (B-13) is treated similarly. For 
this integral, an additional restriction must be included. Within 
these integration limits, Cp = Cp^ sin^O * (1 - 1/4 x^) sin^0, 
and Cp -► 0 as X -► 2, which occurs for angles such that cos 0 -»■ 0. 
This restriction is analogous to the one for normal impingement 
(i.e., r < R < 2 r.) and is caused by a breakdown in the model 
near the periphery of the emergent jet. The numerical result is 


3.T/2 2R^[1 + (1/2 cos 6)] 


-/ / 

tt/2 R. 

3 


V^r dr d0 = ~ [0.172 - 1.576(4 »tt^)^] 


(B-18) 


The average of V^/V^ is the sum of equations (B-17) and (B-18) . 
Next, the average value of AV/V . is computed, again by similar 
numerical integration. The result is 



j/av 


7T^ TT 


-0.281 + 




3 \ 4tt^ 


(B-19) 


Altogether, the average value of V/Vj is the sum of equations 
(B-17), (B-18), (B-19) and one-half of (B-12). Thus, 


^ jo.81 .f - 6.79(*^^)^ - 2^2 to ' 


4tt + 3 ^ 
3 2 \ 47T 


I 


+ (4>tt )^[2 + 2(4»7T.)^+7r-1 + 2{4.TT )^ri+ ((})7rj^l2.n - p- : 


- ((J>7Tj)j^2(<j)'IT^)^ + TT^ \^ + ((J)TT^)‘ £n 


■« 7T^ 
2 . 2 


TT^ + 2(j)7[^ (^1+ (<{)7T~) ^- 


7T2 + 1 


(B-20) 

Because of the many approximations made in the numerical integra- 
tion, equation (B-20) does not represent the "exact" V /V . as 

av 2 

closely as equation (B-12) does for normal impingement. 
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APPENDIX C. POROUS MEDIUM MODEL 


A two dimensional flow model is examined for simplicity. 

The pressure distribution on the surface z = 0 is assumed to be 

w = i pvj [i - i ( 5^)='] ,C-1) 

where is the jet half-width. The coordinate system is shown 
in the sketch. Equation (C-1) is not an accurate representation 



of the impingement pressure for a plane jet but this is of no 
consequence since only a qualitative model is desired to examine 
relative effects. The flow law for porous media is 



where and can be 
resistances across the 
screen. Conservation 



different to account 
screen thickness and 
of mass, 


(C-2) 

for different relative 
in the plane of the 


9v^/9x + 9v^/9z = 0, gives 


9 ^P 


9 x‘ 




9 z' 


4 '3 


+ 


0 


(C-3) 


The boundary conditions are 


V 

P 




2 = 0, X > 0 


(C-4) 


P » P 

o 


z = -b # X > 0 (C-5) 


iz 

9 X 


0 


X = 0, -b <2 <0 (C-6) 


The fourth condition, namely P P as x -b < z < 0 , is 

o “ — 

replaced for convenience by 


P = P X = 2W., -b < 2 < 0 (C-7) 

o 3 - - 


This approximation will overestimate the flow within the screen 
(and underestimate the through flow) since in reality P > P 

o 

at X = 2W . . 

3 

The solution of equation (C-3) which satisfied equations 
(C-5 , (C-6) , and (C-7) is 


00 

= P + / A cosl 

O n L 

n = l 


( 2n - 1 ) TTx 

4W . 

3 


sinh 


(2n - l)7r (z - b) 


L 4W. (K./K- ) 
3 2 1 


1/2 


(C-8) 


The constants A must be determined from boundary condition 
n 

equation (C-4), by expanding [1 - (x/Wj)^/4] into a Fourier 
series of cos [ ( 2n-l ) ttx/ 4W^ ] terms. The result is 




32(-l) 


n+1 


7r (2n-l) sinh 


V 


• r (2n-l)7Tb 1 } 
[4W.(VKj)H *) 


(C-9) 


A value of b/2Wj corresponding to typical tests of reference 4 
is 0.015 cm/0.64 cm = 0.023. Several cases will be examined. 


4 4 


0. This corresponds to straight-through flow for 



which 9 P/3 z = (P - P )/b. Thus 



(C-10) 


(C-11) 


= 1. This corresponds to equal flow resistances in 
the X and z directions. Numerical results are 



(C-12) 


(C-13) 


There is less than 0.07% difference between this flow and 
the one for straight-through flow. 

K^/K^ = 1000. The results are 


0.367 pv. 

3 


-0.0013 pV^ 


(C-14) 


V 

av 



(0.367 + 0.0013 + 



(C-15) 


The through-flow is 80 percent of that for when there is no flow 
in the plane of the screen. 

K^/K^ = 10000. The result^ are 


* 0.0279 PV^ 


A 2 = “0.75 X 10 "^ pV^ 


V =1.39 
av 



(C-16) 


(C-17) 


The through-flow is only 19.2 percent of the flow when no in-plane 
flow exists. Thus, a reduction in flow by a factor of 2, which is 
needed to correlate the data of reference 4, requires a value of 
K^/K^ greater than 1000 but less than 10,000. 


APPENDIX D 


SYMBOLS 


a 

A 

b 

B 

C 

P 

D 

AP 

Q 


r ,0 ,2 


R . 

3 

V 


V 

av 


V , 
D 

V 

s 


X 

a 

6 

e 

w 


TT 


2 


7r 


3 


V' 


surface area to unit volume ratio of screen , cm 

2 2 -1 

screen geometrical constant, aa bQ/E , cm 
thickness of screen, cm 

2 

screen geometrical constant, 6bQ/E D 
pressure coefficient 

screen pore diameter, cm 

2 

pressure drop across screen, N/cm 
tortuosity factor 

coordinate system centered at jet stagnation 
point; r,0 in plane of screen, z positive upward, 
impinging jet radius, cm 
through-flow velocity, cm/sec 
average through-flow velocity, cm/sec 

impinging jet velocity, cm/sec 

boundary layer slip velocity, cm/sec 


viscous resistan 
inertial resista 
screen-volume vo 
liquid viscosity 
blockage coeffic 
dimensionless pa 

dimensionless pa 

dimensionless pa 

dimensionless st 


ce coefficient 
nee coefficient 
id fraction 
, g/cm-sec 
ien t 

rameter, Ay/pv . 

3 

rameter, tt^/B 
rameter, 0.074 
ream function 


(pb V^/yR^ 


1 / 2 , 


1 
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